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Abstract 

We discussed methods used by the PHENIX to constrain the flow background in the two particle jet correlation. 
Both the background level and elliptic flow can be reliably decomposed from the jet contribution. We also studied the 
non-flow contribution to the reaction plane elliptic flow due to dijets. We found the jet bias is negligible in PHENIX, 
when the reaction plane is measured at Beam Beam Counter acceptance (3 < I77I < 4). 



1. Introduction 

Two particle azimuth correlation is a useful tool 
to study the strongly interacting partonic matter 
believed to has been created at RHIC Au+Au col- 
lisions. The typical correlation function can be de- 
composed into a jet part J{A(j)) and an underlying 
event part ^ that is modulated by the elliptic flow [1] : 



C(A(/)) = J(A(/)) + ^(1 + 2?;*w^ cos 2A(/i) 



(1) 



In the case that dijet is not modified by the medium, 
J(A0) behaves just like in p+p collisions. The RHIC 
data has shown a rich modification pattern, which 
is dependent on the pt of both particles and can 
be characterized into four qualitatively distinct re- 
gions. 1) A broadening jet shape at the away side 
and enhancement of the jet multiplicity at both the 
near and away side at low px [2]; 2) A flat or pos- 
sible volcano-shaped away side jet pairs at interme- 
diate Pt [1]; 3) A seemly complete disappearance 
of away side jet or equivalently a flat away side ex- 
tended to the near side at moderately highpT [3];4) 
Reappearance of the away side jet peak at very high 
Pt [4] . Several models were proposed to explain var- 
ious aspects of the modification pattern, but so far 
no model can consistently describe all four regions 
on a quantitative level. 



The medium effects on the dijets can be quantified 
by /aa, which is the ratio of the jet yield per trigger 
in A+A collisions to that in p+p collisions. Since the 
modification is always on the jet pairs, the following 
relation holds between the per-trigger yield using 
the high Pt particles as triggers (type a) and that 
using low Pt particles as triggers (type b): 



-"•aa-'aa — ^aa-'aa — 



JetPairsAA 
iVcoii X JetPairs^ 



(2) 



where i?AA is the single particle suppression fac- 
tor, JetPairsAA and JetPairSpp represent the aver- 
age number of jet pairs in one A+A collision and one 
p+p collision, respectively. One would expect /aa < 
i?AA for hadron-hadron correlation due to trigger 
surface bias. In contrast, the medium is transparent 
to the leading photons in direct photon-jet correla- 
tions. One expect that the away side jet behave ex- 
actly as the single jet suppression: /aa — Raa- 

The level of accuracy in extracting the jet sig- 
nal depends on how well one determines the elliptic 
fiow of both particles, w|, and combinatoric back- 
ground level ^. PHENIX measures the elliptic flow 
through the reaction plane (RP) method, where the 
event plane (EP) is determined by the BBC at for- 
ward region (3 < I77I < 4). The systematic error on 
the V2 is dominated by the EP resolution, 5V2/V2 = 
SV2/V2 — Ercso- £, represents the ratio of the combi- 
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natoric background in same event to that in mixed 
event: {ntUa) = C (^t) (n-a)- It is typically slightly 
above 1 due to finite multiplicity fluctuation in a 
typical centrality bin. The uncertainty on the corre- 
lation function can be expressed as: 

5C^2 {(5vlv^ + C<5u2 ^2 + S^v^v!^) cos 2A0 + 

w (2e,.cso^ + <5C) 2vlv^ cos 2A0 + 6^ (3) 

Given the typical values of various factors: ^ « 1, 
5^ < 0.01, 2vlv^ < 0.1, and e,.oso ~ 0.05 - 0.1, the 
above formula can be simplified as 



5C w 4ei.csoi'2^2 cos 2A(j} + 5S, 



2. Constrain Background Level ^ 



(4) 



PHENIX has previously used two methods in de- 
termining the ^ : the absolute normalization ( ABN) 
method [5] and Zero Yield At Minimum (ZYAM) 
method [1,6]. The ZYAM method assumes there is 
a point in A(j) where jet yield is zero, i.e the back- 
ground term "kisses" the correlation function. The ^ 
value from ZYAM method is a upper limit since the 
background can't be higher than the total in Eq.l. 
The ABN method calculate ^ directly by assuming a 
certain relation between multiplicity and centrality. 

Here we discuss a third method in constraining ^, 
based on the combined information from opposite- 
sign charged pairs (OSC) and same-sign charged 
pairs (SSC). Fig. la shows the azimuthal distribu- 
tions for OSC and SSC for top 0-5% Au-hAu colli- 
sions for 2.5 < PT.trig < 4 GeV/c and 1 < pr.asso < 
2 GeV/c. The OSC correlation at the near side is 
larger than that for SSC, while the charge combina- 
tion has no effect on the back-to-back correlations. 
The difference of the jet strength at the near side 
is a consequence of the charge ordering effect in the 
jet fragmentation process, which leads to more OSC 
pairs than SSC pairs. One would not expect any 
charge correlation between the pair at the Acf) « 
IT, since they come from different jets. The ZYAM 
method clearly finds different ^ values between the 
two although their true values should be identical. 
Fig. lb shows the jet signal estimated independently 
for OSC and SSC using the ZYAM approach. The 
estimated away side jet yield for OSC and SSC pairs 
are clearly different due to the difference in ^ from 
ZYAM. 

In both cases, the ZYAM minimum is at around 
the same A(/)o « 0.8. The difference at that point 
indicated different amount of near side jet contribu- 
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Fig. 1. a),c) The correlation function in 0-5% centrality bin 
for opposite-sign pairs (open circles) and same-sign charged 
pairs (closed circles). The lines indicate the background level 
determined via ZYAM method (a)) and charge-dependent 
method (c)). b),d) The background subtracted distribution 
for using ZYAM (b)) and charge-dependent method (d)). 

tion. Given that the near side jet width are the same 
between SSC and OSC pairs [7] , jet contribution can 
be written explicitly for the two cases as. 



jS = 4(A(/.)-K Ja(A0) 
J^(A0) + J^(A0) 



(5) 



J' 



Define A^ = (^o 



^o4(A0) + Ja(A0) 



i.e. the difference 



in the ^ value between OSC and SSC pairs, we got 



4(A0o) = 



Ao-1 



AC 



^^(0)/4(0) 



^A^ 



jO(A0o) = Ae 



4(0) 



4(0) 



4(0) 



4(0) 



4(0) 



4(0) 



(6) 



(7) 



The near side jet contribution at the A0o can be 
calculated from the near side jet peak value (at A<j> = 
0), which is given by ZYAM method. Since the jet 
fraction is large at the peak region, fractional error 
of the peak values due to uncertain on ^ from ZYAM 
is small. We include it in final systematic error. The 
results of this procedure is shown in Fig. Id. The 
subtracted distribution are identical on the away 
side but shifted above as expected. 

The charged dependence method is useful when 
the statistic is good, such that ^ is well constrained 
at the ZYAM minimum and JN{A(f)o) can be relia- 
bility estimated. We recognize that there could be 
contributions from away side jet at ZYAM minimum 
also, Jj^(A(f)Q), which is inaccessible in current ap- 
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Table 1 

5 values calculated for 2.5 < Pt, trig < 4.0 GeV/c and 1.0 < 
Pt asso < 2.0 GeV/c bin compared between three different 
normalization method for several centrality selections. 



Contrality 


ZYAM 


absolute 


charge dependent 


0-5% 
20-30% 
50-60% 


1.0018 ± 0.0004 
1.015 ±0.0015 
1.076 ± 0.009 


1.0023 ±0.0002 
1.012 ±0.003 
1.07 ±0.02 


0.998 ± 0.002 
1.004 ±0.006 
1.054 ±0.009 



proach. Thus the ^ value obtained in this method 
could still be too big. Table. 1 summarize the ^ val- 
ues obtained from the three methods. 

3. Constrain v-^ from Reaction-Plane 
Dependent Correlation 

Let's define ct„ = (cosn(\E'EP — 4'rp)} as EP res- 
olution for nth order harmonics using the elliptic 
fiow RP. According to [8] , the pair distribution when 
the triggers are selected in a limited angular bite of 

± c with respect to reaction plane is 



Jc rf,,(A0) + 2w"- cos2A(/)) 
a 



(8) 



Sill 2c 

a = 1 + 2vl cos 2(f>s—^ — (72 

f sin 2c f , , sin 4c 

b — + cos 205 (T2 + V2 cos 40s 0-4 

2c 4c 



Fig. 2a shows the CFs in 0-5% central Au+Au col- 
lisions for 6 angular bins in 15° steps. Fig. 2b shows 
the jet yield after subtracting the flow terms calcu- 
lated according to Eq.8. Although the CFs change 
dramatically from in plane to out of plane, the cal- 
culated flow term tracks the true flow background 
nicely. Given the small eccentricity in 0-5%, we can 
safely assume that the jet yield, Jc,4>s does not de- 
pends on the trigger direction. In this case, the small 
difference between the jet functions in Fig. 2b can be 
used to further constrain the V2- 

4. Non-flow Effects from Jets 

Reliable extraction of the jet signal requires accu- 
rate determination of and V2- To this end, con- 
tributions from non-flow correlations that lead to 
azimuth correlations not related to the true RP di- 
rection, need to be studied. The non-flow correla- 
tions include transverse momentum conservation ef- 
fects, resonance decays, HBT correlations that are 
important at low px [9,10] and jet correlations that 
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Fig. 2. a) Correlation function for various 6 trigger direction 
bin and the trigger integrated bin (the center curve, b) The 
background subtracted per-trigger yields, the insert figure 
shows the 6 trigger bins. 

are important at high px [11]. The non-flow correla- 
tions affect the V2 values either by changing the EP 
resolution or cause fake V2 by biasing the EP direc- 
tion towards the now-ffow particles. The former is 
small if the non-flow particle multiplicity is small. 
The later can be suppressed in PHENIX if the cor- 
relations are limited in a narrow Arj window, such 
as the resonance and intra-jet correlations, thanks 
to the large rj separation between BBC and central 
arm. However, the inter-jet correlations can still bias 
the BBC reaction plane determination, due to their 
broad distribution in Arj. 
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Fig. 3. The effect of the dijet on reaction plane resolution 
as function of centrality in different rapidity windows. The 
embedded jet is at mid-rapidity. 

We study the biases due to dijet correlations 
by embedding dijet pairs into background events 
with realistic flow modulation. The background 
Au+Au events are simulated with HIJING, which 
was checked to reproduce the charged hadron mul- 
tiplicity in rj from PHOBOS. Elliptic flow is imple- 
mented by applying a track by track weight for each 
HIJING event: 



w{b,pT,ri) = l + 2v2{b,pT,r]) cos2{ 



where the ^ is the direction of the impact parameter 
b. The centrality, pt dependence of the V2 is tuned 
according to the PHENIX measurement [12]. We 
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used a common 77 dependence from PHOBOS [13] 
for all centrality selections. The dijet pairs are gen- 
erated from PYTHIA event generator, requiring 
a leading particle above 6 GeV/c at mid-rapidity 
(|r/| < 0.35). This corresponds to typical energy of 
6/(z) « 10 GeV/c for the dijets. 

The effects of the dijets are evaluated by compar- 
ing the event plane before and after the embedding. 
Since dijet pairs are random with respect to the RP, 
the EP resolution for combined event is worse as 
shown in Fig. 3. The difference in EP resolution is 
small except in peripheral bins where the event mul- 
tiplicity is small. And the difference become negligi- 
ble in more forward rj window. On the other hand, 
dijet tends to pull the event plane towards the di- 
jet direction, resulting in a fake V2 for the leading 
hadrons as shown in Fig. 4, where the relative az- 
imuth distribution between the leading hadrons and 
EP from either the HIJING event or the combined 
event are plotted. The dijets are clearly correlated 
with the EP determined from the combined event, 
thus have a fake V2 ■ 




"o 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 
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Fig. 4. The distribution of the leading particle from the 
dijets relative to the event plane calculated from HIJING 
only (loft) and event plane from the embedded event (right). 

The size of the fake V2 depends on the rapidity 
gap between the triggering hadron and the subevent 
used to determine the event plane. Due to the away 
side jet swing, one would expect the this bias persists 
to large rj region. Fig. 5 shows the centrality depen- 
dence of fake V2 for various rapidity window. This 
fake V2 is the raw signal extracted from fit in Fig. 4 
without divided by EP resolution. The fake V2 de- 
creases as the subevent moves towards large 77. When 
the subevent is 3 < \ri\ < 4 (BBC acceptance), the 
fake V2 becomes almost negligible. 

So far, we have assumed normal dijet in our study. 
RHIC data indicates that there is a broadening of 
the away side jet and an increase of jet multiplicity 
at low pt by about factor of 2 [2,18]. To account for 
that, we increase the PYTHIA dijet multiplicity by 
factor of 2 and redo the study of the rapidity depen- 
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Fig. 5. The fake V2 of the leading particle as function of 
centrality using the EP determined in four rj windows. The 
embedded jet is at mid-rapidity. 

dence of the fake V2 ■ We did not take into account the 
broadening of the near side jet in 77, since it is dom- 
inated by the away side jet swing. Fig. 6 shows the 
fake V2 corrected for the EP resolution as function 
of rapidity window with enhanced dijet and normal 
dijet. The elliptic flow is higher when the jet mul- 
tiplicity is doubled, but it is still negligible in the 
BBC acceptance. 
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Fig. 6. The EP resolution corrected fake V2 as function of 
^low < I*?! < 2.8 for normal dijet and enhanced dijet. 

The V2 measurements can also be affected by the 
event by event fluctuation of V2 [15,16,13,17], which 
is a direct consequence to the event by event fluctu- 
ation of the collision geometry. However, since the 
two particle correlation method is used in this anal- 
ysis, the event by event fluctuations would also con- 
tribute to the V2 correlations. In this sense the W2{2} 
should be the one to use in Eq.l. Since i'2{2} was 
shown to be consistent with the BBC RP V2 [14] up 
to 4 GeV/c in pT, it is reasonable to using BBC RP 
V2 in Eq.l. 
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